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A microscopic analysis of the superconducting quantum critical point realized via a pair-breaking 
quantum phase transition is presented. Finite temperature crossovers are derived for the electrical 
conductivity, which is a key probe of superconducting fluctuations. By using the diagrammatic 
formalism for disordered systems, we are able to incorporate the interplay between fluctuating 
Cooper pairs and electrons, that is outside the scope of a time-dependent Ginzburg Landau or 
effective bosonic action formalism. It is essential to go beyond the standard approximation in 
order to capture the zero temperature correction which results purely from the (dynamic) quantum 
fluctuations and dictates the behavior of the conductivity in an entire low temperature quantum 
regime. All dynamic contributions are of the same order and conspire to add up to a negative total, 
thereby inhibiting the conductivity as a result of superconducting fluctuations. On the contrary, the 
classical and the intermediate regimes are dominated by the positive bosonic channel. Our theory 
is applicable in one, two and three dimensions and is relevant for experiments on superconducting 
nanowires, doubly-connected cylinders, thin films and bulk in the presence of magnetic impurities, 
magnetic field or other pair-breakers. A window of non-monotonic behavior is predicted to exist as 
either the temperature or the pair-breaking parameter is swept. 
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I. INTRODUCTION 

The study of quantum phase transitionsiiS. has consis- 
tently been one of the frontier fields in condensed mat- 
ter physics, for the last decades. The steady growth in 
the possibility of simultaneously accessing lower temper- 
atures and higher values of tuning parameters such as 
pressure and magnetic field, together with the possibility 
of operating old and new experimental techniques un- 
der such extreme conditions, has provided the necessary 
thrust. Discovery of new materials and physical systems 
has been another key factor. Not only has it maintained 
the freshness and novelty of the field but also helped to 
identify the universal features. Complex materials usu- 
ally go hand in hand with complex phase diagrams with 
many competing or coexisting orders, given the multitude 
of energy scales that they have. Although this accounts 
for a rich body of physics, it also makes it harder to un- 
mask the universal features and to systematically explore 
the neighborhood of the quantum phase transitions in the 
system. 

Technological advances in relatively recent years have 
made it possible to make precisely designed and control- 
lable bulk, mesoscopic and nano systems, consisting of 
cold atoms and quantum dots, for example. These have 
not only opened new avenues for studying quantum phase 
transitions but also served as model systems where theo- 
retical predictions can be verified with the help of tunable 
parameters. This in turn can act as a starting point for 
understanding complex materials which are not as easy 
to control. 

By now there is a large body of theoretical work de- 



voted to studying quantum phase transitions in a whole 
range of systems ranging from heavy fermion compounds, 
high temperature superconductors, manganites and or- 
ganic materials, to quantum dot and cold atomic sys- 
tems. Transitions between and out of different correlated 
states of matter ranging from magnetic, superconduct- 
ing, and charge-ordered states to more exotic, fractional 
and topological states, have been the subject of research. 
In spite of active interest, many questions still remain 
unanswered. 

To-date superconductivity remains one of the most 
striking examples of emergent many-body states and 
quantum phase transitions involving a superconducting 
state are ubiquitous in a whole variety of materials. In 
many a cases, it appears as one of the multiple phases in 
a complex phase diagram, with its mechanism not always 
clearly understood. However there is another set of ma- 
terials consisting of single elements or simple compounds 
with relatively simple phase diagrams which too display 
superconductivity. 

The BCS theory that was put forward exactly fifty 
years ago was able to explain the basic mechanism of 
superconductivity which in turn was discovered for the 
first time almost another fifty years earlier. This the- 
ory is one of the crowning glories of condensed matter 
physics and a prototype of an ideal many body theory, 
which is successful in explaining the experimentally ob- 
served thermodynamic and electrodynamic behavior in 
a large class of materials which is commonly labelled as 
conventional or BCS superconductors^. 

The point of view we want to take in this paper 
is to study the physics of quantum phase transitions 
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in such a BCS superconductor for which the theory 
is well-understood and experimental properties well- 
characterized. The theory could then be cleanly tested 
by taking a conventional superconductor with a simple 
phase diagram. Such a study would exemplify the physics 
of a superconducting quantum phase transition which is 
not only of fundamental and technological interest but 
also a prototype for quantum phase transitions in other 
correlated systems. 

The key characteristic of the BCS theory is the pairing 
of electrons with their time-reversed partners to form a 
condensate of Cooper pairs that superconducts. As the 
temperature is increased more and more quasi-particles 
are formed by exciting electrons out of the condensate by 
breaking the pairs, until superconductivity is destroyed 
at a mean-field transition temperature. The experimen- 
tal transition temperature is well-defined and the transi- 
tion is mean-field like in most conventional bulk super- 
conductors. Superconducting fluctuations however do ex- 
ist even beyond the transition temperature and the fluc- 
tuations become stronger as the dimensionality of the 
sample is reduced. The effect of fluctuating Cooper pairs 
on different physical properties such as diamagnetism, 
specific heat, and conductivity, was actively studied - 
both theoretically and experimentally- about a decade 
after the BCS theory was put forward; see the review 
on fluctuation effects in Ref. [5| which was already writ- 
ten three decades ago. The interest in studying fluctu- 
ation effects saw a revival after the discovery of high- 
temperature superconductivity and corresponding the- 
oretical results can be found in a more recent review 
article^. 

What is the way to start from a superconducting state 
at the absolute zero of the temperature and destroy it 
at a flnite value of some tuning parameter, via a second 
order phase transition? One has to think of a way of 
breaking the pairs without the help of thermal fluctua- 
tions. It turns out that all one has to do is to turn on 
a perturbation that breaks the time-reversal symmetry. 
Pair-breaking perturbations resulting in a suppression 
of the transition temperature have been well-understood 
and it has been known for a long time that supercon- 
ductivity is destroyed even at zero temperature, once the 
pair-breaking parameter reaches a critical strength. We 
recognize such pair-breaking quantum phase transitions 
(with dynamic critical exponent z = 2) out of a super- 
conducting state as an important class of quantum phase 
transitions and make it the subject of our study. 

In spite of the long history of superconducting fluctu- 
ation effects near the classical transition, there has been 
no systematic theoretical analysis of such effects in the 
vicinity of a pair-breaking quantum phase transition. In 
this paper we have used the flnite temperature Matsub- 
ara diagrammatic method to evaluate the fluctuation cor- 
rections to the electrical conductivity which is one of the 
key physical quantities capable of diagnosing and show- 
casing the role of superconducting fluctuations. We have 
incorporated the (dynamic) quantum fluctuations by go- 



ing beyond the standard approximations and mapped out 
the fluctuation regimes near the transition. Such a study 
is timely given the progress in fabricating ultra-narrow 
superconducting nanowires and doubly-connected cylin- 
ders in addition to thin-fllm samples. Not only will it 
allow identifying the universal features near supercon- 
ducting quantum phase transitions in complex materials, 
but also play a crucial role in enhancing our understand- 
ing of mesoscopic superconductivity. The microscopic 
approach we use allows us to treat the fluctuation cor- 
rections involving the interaction between the electrons 
and the fluctuating Cooper pairs and identify the regimes 
in which they dominate. Our analysis can then serve as a 
guideline for the construction and validation of an effec- 
tive bosonic theory or a time-dependent Ginzburg Lan- 
dau formalism which does not incorporate the electrons. 
Such a bosonic theory can then be used for further anal- 
ysis of the quantum critical point especially in one and 
two dimensions, given that the upper critical dimension 
in our case is two. 

The outline of the paper is as follows: In Sec. II we 
lay the foundation for the remainder of the paper, by 
clearly defining the meaning of a pair-breaking perturba- 
tion and introducing the class of pair-breaking quantum 
phase transitions that are of interest to us. We present 
the framework based on Usadel equations as a systematic 
and general method for obtaining the expression for the 
pair-breaking parameter in a given situation and illus- 
trate it by considering some examples which are simple 
yet relevant for our analysis. 

In Sec. Ill we present the calculation of the fluctuation 
corrections to the normal state conductivity in the vicin- 
ity of the pair-breaking phase transition. We introduce 
the key building blocks of the temperature diagrammatic 
perturbation theory that we need, show how they are 
modified in the presence of a pair-breaking perturbation, 
and present their limiting forms. A fairly detailed ac- 
count of the actual evaluation of diagrams corresponding 
to the fluctuation corrections is then given, focusing on 
the careful considerations required in order to incorpo- 
rate the quantum (dynamic) fluctuations correctly. The 
results of this section are applicable near the entire phase 
transition line, starting from the classical flnite temper- 
ature transition in the absence of pair-breaking pertur- 
bation to the quantum phase transition driven by tuning 
the pair-breaking parameter at zero temperature. 

The evaluation of these general expressions in the 
vicinity of the quantum phase transition is done in Sec. 
IV. Once the dominant corrections are identified, we 
present the different fluctuation regimes that come out 
of our analysis. There are three regimes -quantum, in- 
termediate and classical- and the behavior of the conduc- 
tivity depends on the path of approach to the quantum 
phase transition. We demonstrate how these results vary 
based on the effective dimensionality of the problem and 
provide predictions of our theory that should be appli- 
cable to experiments on nanowires or doubly-connected 
cylinders, thin fllms, and bulk systems. The result we ob- 
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tain by evaluating the "Aslamazov-Larkin" correction in 
the vicinity of the classical transition is given as a bench- 
mark to compare with the well-established results in the 
literature and also with its behavior near the quantum 
phase transition. 

In the final section we place our work in a bigger pic- 
ture by discussing related theoretical and experimental 
work. Different theoretical approaches for the same prob- 
lem as well as slightly different physical configurations 
studied using a formalism similar to ours, are both in- 
cluded. An attempt to interpret all the experiments on 
superconducting quantum phase transition in thin films 
would clearly be outside the scope of this paper. We have 
hence focussed mainly on those that have analyzed the 
data in terms of quantum corrections, including the fluc- 
tuation corrections. Relevant experiments on supercon- 
ducting nanowires and doubly-connected cylinders are 
relatively scant so far. However we discuss the current 
status to support our belief that given the technologi- 
cal advances in the recent years, the predictions of our 
theory should not only be accessible but also important 
from the point of view of using ultra narrow wires in 
superconducting electronic circuits. 



II. PAIR-BREAKING PARAMETER (q) 

A. Definition and physical meaning 

A pair-breaking perturbation is any perturbation that 
breaks the time-reversal degeneracy of a superconduct- 
ing paired state. Anderson's theorem^ asserts that in 
the absence of such a perturbation, the superconducting 
critical temperature and the BCS density of states re- 
mains the same even after alloying the superconductor 
with impurities. However if the impurities are magnetic, 
Abrikosov and Gorko\«^ found that the Tc is suppressed 
and the density of states is modified as well (giving a 
gapless regime) . They parametrized the strength of the 
pair-breaking perturbation by a pair-breaking parameter 
a and obtained 
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where ^ is the digamma function and Tcq = Tc{a — 0). 
The parameter a was shown to be inverse of the spin-flip 
scattering time which is proportional to the density of 
magnetic impurities. 

Following this classic work, it was recognized that the 
above equation for Tc suppression (as well as gapless su- 
perconductivity) can be transcribed for a whole class of 
pair-breaking perturbations for which the transition to 
the normal state in the presence of a is of second order, 
once the appropriate a is used for each case^i^. What is 
essential is the presence of a rapid scattering mechanism 
that modulates over time the pair-breaking perturbation 
seen by a given Cooper pair of electrons, assuring an er- 
godic behavior of electrons. Then a can be interpreted 



as the depairing energy (energy splitting) of a pair of 
time-reversed electrons, averaged over a time interval tk 
it takes for their relative phase to be randomized by the 
perturbation; one thus has 2a = H/tk- 

This generalization of the concept of pair-breaking of- 
fers the possibility of defining a pair-breaking parameter 
not only for bulk systems but also for mesoscopic and 
non-homogeneous systems and for situations in which 
multiple pair-breaking mechanisms might be operative. 
However, the derivation of a might not be straightfor- 
ward in such cases. Below we will describe the Usadel 
equation formalism as a general method to derive a for 
a given configuration. 



B. Derivation using Usadel equations 

In this paper we are interested in focussing on 
dirty superconductors for which the Usadel equation 
formalisrri^ii^ is well-suited. Writing the Heisenberg 
equation of motion starting from the BCS Hamiltonian, 
one gets the microscopic Gorkov equations for the normal 
and anomalous Green functions. Using the fact that the 
characteristic length scale for the normal state is smaller 
than the length scale for the superconducting order pa- 
rameter variation, one can make the quasiclassical ap- 
proximation 
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(which is quite accurate for most classic low-temperature 
superconductors and less so for high temperature super- 
conductors) to exclude the fast oscillations of the Green 
functions associated with variations of the relative coor- 
dinate on a scale h/pp and rewrite the Gorkov equations 
in terms of the quasiclassical Green functions with only a 
slow dependence on the centre-of-mass coordinate vary- 
ing on the scale ■ The Eilenberger equations so derived 
can be further simplified in the dirty limit 
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(where I is the mean free path and t is the mean free 
impurity scattering time) in which the strong scattering 
by impurities produces averaging over momentum direc- 
tions, to obtain the Usadel equations^ 

-iD[g{V - —Aff - /V^g] = 2Ag - 2zc^„/, (4) 
c 

-iD[g{W + —Afp ~ fW'g] = 2A*g - 2*c.„/t (5) 



where D = VpT/3 is the three-dimensional diffusion con- 
stant. Note that g and /, the quasiclassical Green func- 
tions averaged over momentum directions, get expressed 
as functionals of the fluctuating order parameter field 
A{x,t). Below we outline the precise recipe for deriving 
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the pair-breaking parameter, starting from these equa- 
tions. 

To identify the pair-breaking parameter it is enough to 
consider the Usadel equations for the / function to the 
lowest order in A, so as to obtain 

-iDgoiV -—Aff = -2iLuf + 2Ago (6) 
c 

with go — sign a; . One could interpret the depairing pa- 
rameter as the lowest eigenvalue of the operator (in the 
transverse direction) 



with 



af 



(7) 



obtained by appropriately choosing the gauge and the 
boundary conditions (to ensure absence of current per- 
pendicular to a wire or film surface, for example). Al- 
ternatively one can solve Eq. ([6]) for / and read out the 
expression for a from it. In the remaining part of this 
subsection we illustrate the procedure for some examples 
of interest. The formalism is general enough and can be 
appropriately adapted to complicated situations involv- 
ing the simultaneous occurrence of multiple pair-breaking 
perturbations. 

Although the pair-breaking parameter for the classic 
example of magnetic impurities can be derived within 
the Eilenberger-Usadel formalism, we refer the reader to 
the original paper— and proceed to discuss some other 
specific examples. 
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identified as the pair-breaking parameter. 



2. Nanowire 

Consider a wire of radius R such that the diameter 
is smaller than ^ and A and the effective dimensionality 
of the problem is d = 1 as far as the superconducting 
fluctuations are concerned. To obtain the expression for 
the pair-breaking parameter coming from the orbital ef- 
fect of a magnetic field H applied parallel to the wire, 
we use the cylindrical coordinates and choose the gauge 
such that A^f, = Hp/2, Ap — = 0. For the sake of 
illustration, this time we start with Eq. ([7]) which takes 
the form 
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and on integrating in the transverse direction (i.e. over 
p) gives 

R ^ Ft 

^f[pdp(^\ ^affpdp (13) 



1. Thin Film 

Consider a thin film with thickness s smaller than the 
superconducting coherence length ^ such that the super- 
conducting fluctuations are effectively two-dimensional 
(d = 2). Let us focus on the orbital pair-breaking effect 
of a magnetic field H applied parallel to the film. We 
consider s to be smaller than the penetration depth A so 
that one can assume the field to be uniform inside the 
sample. If the film is placed in the xy— plane and H is 
parallel to the y— axis, we choose the gauge such that 
Ax = Hz, Ay = Az =0, where z is measured from the 
mid-plane of the film. Now Eq. ([6]) takes the form 
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f = -i\Lu\f + A 



(8) 



which on integrating over the transverse direction (i.e. 
over z) leads to 
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on using the appropriate boundary conditions that re- 
quire the transverse derivative of / to vanish on the sur- 
face of the film. We thus have 
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using the boundary conditions that require df /dp to van- 
ish at p = R. The expression for the pair-breaking pa- 
rameter 
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is then immediately evident. 

If one considers a magnetic field applied perpendicular 
to the wire, the pair-breaking parameter is given by 
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instead and the calculation follows on the lines similar to 
that for a field applied parallel to a film. 



3. Doubly-connected cylinder 

Consider a doubly-connected (hollow) cylinder with in- 
ner radius ri and outer radius r2 such that the wall of the 
cylinder is thinner than ^ and A and the effective dimen- 
sionality is d = 1 as for the case of a nanowire. Due to 
the single-valuedness of A and / their 0— dependence is 
given by e'""^ where n is an arbitrary integer. Choosing 
the cylindrical gauge 
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as we did also for the nanowire, Eq. ^ becomes 



Pair-breaking phase transition 
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Integrating over the radial direction (from ri to r2), we 
obtain the expression 
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c 



,ln(r2/ri) 



(18) 

where n is an arbitrary integer (note that for ri = 
we correctly recover the result for the nanowire). For a 
thin cylinder (ri « r2 ~ r) the pair-breaking parameter 
reduces to 



a = {D/2r'^){<S>/^o-nf 



(19) 



where $ is the flux enclosed by the cylinder, thereby ren- 
dering the classic Little-Park oscillations'^ of Tc as can 
be seen from Eq. ([T]). Interestingly, for a cylinder with 
small enough radius, r < Tc = \/ Dj/AttTco, it is possible 
to push the Tc down to zero at magnetic fields corre- 
sponding to half- integer fluxes ^ = ^o{l/2 + n). 



1 r doubly-connected 
cylinder 




Superconding phase Normal phase 



0.4 0.6 



QCP 



FIG. 1: Phase diagram showing the pair-breaking transi- 
tion from superconducting to normal state with the bound- 
ary given by Eq. pO)) . A superconducting quantum critical 
point is seen to be realized when the pair-breaking parameter 
reaches a critical strength aco/Tco = 0.889. Different systems 
of interest are also illustrated. 



C. Pair- breaking phase transition 

The boundary between the superconducting and nor- 
mal state in the a — T plane is given by 
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At a given pair-breaking strength a, superconductivity 
is destroyed at T = Tf,(a) and at a given temperature 
T, at a = ac(T), obtained by solving Eq. ^ for T 
and for a, respectively. In the absence of a pair-breaking 
perturbation (a = 0) the system undergoes the classi- 
cal transition at Tc(0) = Tco- In the neighborhood of 
this classical transition, for a <^ Tco we can define the 
quantity 



OT[a,T) = 

Tc{a) 



(21) 



that measures the relative distance from the critical tem- 
perature Tc{a). 

On the other hand, if the pair-breaking effect is suf- 
ficiently strong, superconductivity is destroyed even at 
T = thereby yielding a second order quantum phase 
transition. The critical value of a at zero temperature. 



adT = 0) = 



ttTco 
~2T 



can be obtained by using the first term in the asymp- 
totic form of the digamma function for large arguments 
{iIj{z) = Inz). Here ln7 « 0.577 is the Euler constant 



and one can immediately see that 2acO — l-76Tc — Aq, 
the BCS gap at zero temperature. Expanding the RHS 
of Eq. ((20|) to next order in small T, one finds that the 
transition curve close to aco is given by 



ac{T < ac 



6acO 



acO 



(22) 



In the vicinity of the the pair-breaking quantum phase 
transition we will define the quantity 



Sa{a,T) 



«c(r) 



«c(r) 



(23) 



which can be interpreted as the relative distance from the 
critical pair-breaking strength ac(T) at a given T <^ aco- 
In the following section we will be interested in eval- 
uating the fluctuation corrections to the normal state 
conductivity in the vicinity of a pair-breaking transition. 
The effective quasiclassical approach based on the Us- 
adel equations that we have discussed in the previous 
subsection is quite accurate in the dirty limit and has 
no applicability restrictions in terms of the temperature 
range. To evaluate the fluctuation corrections, one could 
imagine using the functional formalism based on this ap- 
proach. However we have chosen to use the standard 
diagrammatic method instead. 



III. SUPERCONDUCTING FLUCTUATION 
CORRECTIONS TO CONDUCTIVITY 

We will carry out a microscopic calculation within the 
standard framework of temperature diagrammatic tech- 
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nique for a disordered electron system^^'^^ in the diffu- 
sive limit (t~^ » T,a). This technique has been ex- 
tensively used in studying the weak localizatioi*^^ and 
electron-electron interaction^^ corrections to the conduc- 
tivity in low-dimensional systems. In the same way as 
these corrections were studied also using alternative for- 
malisms including the non-linear sigma model^S, it seems 
plausible to have an alternative derivation of the super- 
conducting fluctuations corrections to the conductivity, 
which falls in the same league. Here we will restrict our- 
selves solely to diagrammatic perturbation theory. 



A. Basic ingredients 

Although the framework we use is standard, we will 
briefly discuss all the ingredients we need mainly for two 
reasons. First, we want to precisely demonstrate the way 
in which the presence of a pair-breaking parameter mod- 
ifies these ingredients. Second, we want to catalogue all 
the expressions we need, including their limiting forms in 
the vicinity of the quantum phase transition. 



1. Green function 

As is standard, we assume a random disorder potential 
V{y) drawn from a Gaussian white noise (^-correlated) 
distribution such that (^(r)) = and (F(r)F(r')) = 
{V"^) 5{r — r'). In the diagrams, a dashed line denotes 
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(24) 



where v is the density of states at the Fermi level and , 
as defined earlier, is the frequency of elastic collisions. 
The single electron Green function -denoted by a full 
line in the diagrams of Fig. [21- is given by 



G(tJ„,p) 



1 



i[LOn + 



sign(^„) 
2t 



(25) 



where is the single particle excitation spectrum mea- 
sured from the chemical potential. 



2. Cooperon 

Diffuson and Cooperon are the key correlators, repre- 
sented by a sum over ladder diagrams involving coherent 
scattering by impurities, in the particle-hole and particle- 
particle channel, respectively; the latter being of interest 
to us in the context of superconducting fluctuations. The 
expression for Cooperon -represented by a shaded rect- 
angular block in the diagrams- in the presence of a pair- 
breaking parameter a is given by 



C(u;„,a;m,q) 



2'n:v9{-U!nUJm) 



(26) 



where 9{x) is the Heaviside theta function and 

aq = a + Dq^l2 (27) 

Here a;„ is a fermionic Matsubara frequency and q is 
the momentum in the effective dimension as far as the 
superconducting fluctuations are concerned (note that D 
is the diffusion constant in three dimensions as long as 
the diffusion is still three-dimensional). 

Coherent scattering on the same impurity, by both 
the electrons forming a fluctuating Cooper pair leads to 
renormalization of the vertex part in the particle-particle 
channel, given by 



A(w„,a;m,q) 



2'KVT 



(28) 



In the diagrams shown in Fig. ^ we denote A by a shaded 
triangle. 



3. Fluctuation propagator 

The main building block of the diagrammatic tech- 
nique that encodes the BCS superconducting interac- 
tion is the so-called fluctuation propagator (represented 
by a wavy line in the diagrams). It is the impurity- 
averaged sum over the ladder diagrams corresponding to 
the electron-electron interaction in the Cooper channel. 
The expression in the presence of a, obtained using Eq. 
(|26p in a standard way, is given by 



^"'(l^^.|,q) 
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(29) 

where il^, is a bosonic Matsubara frequency. 

The pole of Eq. (|^^ for g, = traces the boundary 
between the superconducting and normal phases (see Sec. 
II C). In the limit of zero pair-breaking strength (a — 0) 
and near the classical transition, T ~ T^o, one can show 
that the expression for the fluctuation propagator reduces 
to 



K-\\a,\,q)^ST{0,T) + 



Dq 



AttT ^ \2 



(30) 



where -0' (1/2) — 7r^/2 and (5^(0, T) is defined earlier by 
Eq. ^ 

On the other hand, at low temperatures, T <^ ckco, the 
fluctuation propagator can be reduced to a form 



i^-i(|f^,|,q) - In 



\n,\/2 



(31) 



which correctly reproduces the transition curve a — 
ac(T). In certain regimes (see below) it is legitimate to 
expand the logarithm and get an even simpler expression 



K{\n,\,q) 



1 



Saia,T) 
where da{a, T) is given by Eq 



2ao(r) 



(32) 
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"Aslamazov-Larkin" 




"Density of States" 




"Maki-Thompson" 




FIG. 2: Diagrams for the fluctuation conductivity divided 
in three groups. There is one diagram corresponding to the 
positive "Aslamazov-Larkin" correction. Three diagrams cor- 
respond to the negative "density-of-states" corrections, and 
each of them has two possible ways of putting arrows on the 
electron Green functions. There are two diagrams correspond- 
ing to the "Maki-Thompson" interference correction with no 
prescribed sign, the second of which has two ways of putting 
the arrows. Full lines stand for the disorder averaged normal 
state Green's function, wavy lines for the fluctuation propaga- 
tor K, the shaded rectangles for the Cooperon C and shaded 
triangles for the vertex C/2nvT. 



B. Evaluation of the diagrams 

The boundary between the superconducting and nor- 
mal region in the a-T phase diagram is given by Eq. 
(|20p . Even while superconductivity is destroyed, super- 
conducting fluctuations continue to persist in the normal 
state region and modify the normal state conductivity. 
We will evaluate the corrections to the conductivity com- 
ing from superconducting fluctuations, using the stan- 
dard Kubo formalism for linear responseii. The electro- 
magnetic response operator Q{flf_i) is evaluated and the 
external frequency is analytically continued into the up- 
per half-plane of the complex frequency {i^l^ — » fi). The 
fluctuation conductivity can then be obtained using 



Sa^n) = lim 9tl^ 



(33) 



once the appropriate Q{fl^) has been evaluated. 

In what follows we will evaluate Q(r2^) using a stan- 
dard set of diagrams^. As shown in Fig. O these 
diagrams can be divided into three groups based on 
their physical interpretation. Of these, the "Aslamazov- 
Larkin" (AL) type of correction^^ is the most intuitive. It 
is a positive contribution coming from additional charge 
transfer via fluctuating Cooper pairs. The negative "den- 
sity of states" (DoS) correction results from the reduc- 
tion of the normal single-electron density of states at the 
Fermi level after accounting for the electrons participat- 
ing in fluctuating Cooper pairs. The corresponding dia- 



grams have only one electron line affected by the fluctu- 
ation propagator. The third, and the more indirect, cor- 
rection is given by the "Maki-Thompson" {MTy^Ld^ 
diagrams which could be thought to be originating from 
coherent Andrecv scattering off the fluctuating pairs. 
One can see from Fig. [2] that the fluctuation propagator 
brings about interference between the electron lines. 

Let us start by considering evaluation of the 
Aslamazov-Larkin diagram. The effective triangular ver- 
tex on either side (both yield the same final expression) 
is given by 

r(17i^, f2^, q) = eT ^ p^X{-uJn + ^lu, i^n - fi^, q) 

G{-ujn + -p + q)G(cj„ - ^f,,p) 
G(tJ„,p)A(cj„,-cj„ + f7i^,q) (34) 

(having chosen the x direction for concreteness) . The 
presence of Heaviside 9 functions in the expression for 
the vertex renormalizations A, defined by Eqs. (|26p and 
(PS)) . dictates the possible signs and ranges of different 
frequencies. Taking these into account and by making 
simplifications valid in the small q and fli^, <C 
limit of interest to us, we find a rather compact expres- 



(35) 



with 



(36) 



and 



1 a 
^ ' 2 ~ 2t:T 



2 



2ttT 



(37) 

where ipi^) is the digamma function. 

To evaluate the entire AL diagram, let us consider the 
summation over the internal bosonic frequency, 

/(r!^) = Tj2Bini,,nf,fK{\ni,\)K{\ni, - 

(38) 

and write it as contour integration in a standard way 
(note that we have temporarily suppressed the q depen- 
dence in the fluctuation propagator K for the sake of 
compactness). By taking into account the analytic! ty of 
the integrand, the evaluation of / is reduced to an inte- 
gration across two branch-cuts. After combining terms 
and analytically continuing the external frequency to the 
upper half plane (if2^ — s- fl), one flnds that 



with 



in 



la 



^^/v(-zl7i)i^(*r!i) 



(39) 



sinh 



2T 



(40) 
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and 



or o 

h^— dnicoth-^K{~ini-in)K(-ini) 

4:TTi J 2T 



ip{-ini~in,-in) + tij{~ini,-in) (4i) 



Since the contribution to the conductivity (Eq. ([55]) ) 
goes as I/^, we need only consider the terms in / that 
are hnear in Q since we are interested in the 17 ^ limit; 
terms that are zeroth order in fl will be cancelled by 
analogous terms in the remaining diagrams to ensure the 
absence of anomalous diamagnetism in the normal state. 
Carrying out an expansion in 17 we can write 



tpl—ifli, —ifl) 



1 



47rr 

1 in 



1 aq — ifli/2 

2 ^ 



2ttT 



(42) 



2 (47rr)' 



2-kT 



Since la is already linear in 17, in it we need only keep 
the zeroth order term from this expansion. On the other 
hand, using the expansion for both K and ij:, we have 



-17 f 17i 

Ih — / dVLi coth — 

2-K 2T 




We can integrate the first term by parts, combine simi- 
lar terms from la and lb, and after some manipulation 
get a final expression for /. One can then immediately 
write down the Aslamazov-Larkin fluctuation correction 
to conductivity as a sum of two terms: 



6a 



AL 



Sa 



AL 
sh 



6a, 



AL 
cth 



where 



_ D^e^ f dfiq dl7i 
~2^J {27r)<^sinh2^ 

(Im{ifH17i,g)7H17i,q)})' 



+ Im{/f (-zl7i, q) 7(-zl7i, q)'}Im{i^(-zl7i, g)} 



(44) 



(45) 



Going through the derivation, the reader can readily 
convince her/him-self that the contribution (5ct^^ would 
be missed if one were to make the so-called static approx- 
imation in the effective vertex and use r(17ii, — 0, 17^, q) 
in the evaluation of the Aslamazov-Larkin diagram. As 
long as one is interested in obtaining the corrections near 
the classical transition, the approximation is justified and 
6afi^ is indeed the dominant contribution. We will show 
below that in this limit, the result we obtain is in agree- 
ment with the existing literature^. In the a <C T ~ T^q 
limit, we can expand 



7(-il2i,q) 



2ttT 



a. 



where we have used 17i <^ T, in addition. Keeping the 
zeroth order term amounts to making the static approx- 
imation (17i = 0) and Eq. (|45)l reduces to 



= d (2^^ [-2 



d q dfli 2 



g2[Imif(-^17l,g)]2 (49) 



Using Eq. pO|) for the fluctuation propagator in the limit 
under consideration, we have 



lmK{-ini,q) 



17i 



^-'(1/2) {6T + Dq^)^ + nl 



where 



St^6t{0,T) — 



^'(1/2) 



(50) 



(51) 



Further using sinha; — > x for small x, we can write the 
final expression 



6a 



AL_ 
sh 



IGTD^e^ f d'^qdn^ 



nd 

rp2 



{2t:Y 



d'^q 



{6T + Dq'Y + ni 

„2 



d(27r)" 



2T 



(52) 



where we have redefined 



and 



6a, 



AL D^e^i f d'^qdVl^ 



cth 



with 



7(-ilii,q) 



{2TT)'i 

g-»»i/2 
2nT 



coth 



^qlK'i-^n,,q) 

fl aq-^n,/2 
^ \ 2 2ttT 



2t:T^ \2 



il7i/2 



2ttT 



(46) 



(47) 



St — 



4:6t{0,T) _ 6j 



2T 



(53) 



to put the integral in a form similar to what we will find 
in one of the regimes near the quantum phase transition. 
The presentation of the results for different dimensions 
will accordingly be deferred to Sec. IV in order to facili- 
tate the comparison. 

Now we move on to the other two corrections given 
by the "density-of-states" and the "Maki-Thompson" di- 
agrams as shown in Fig. [21 The calculation follows 
on lines similar to what we have outlined for the case 



9 



of Aslamazov-Larkin diagram in detail above. We find 
that the density-of-state fluctuation correction can be ex- 
pressed as 



(54) 



where 



DoS _ ^^e^ f d'^qdmdn uj iK{-in,q) 



AwT J {2t:Y 
-Re ^ 



2T sinh^JI 



Re- 



and 



with 



5ogf = -A+-B 



De^ f df^qdujdVL , LU , 

A — / — ; — -1 — tanh — coth — 

27r / 2tt¥ 2T 2T 



1 



:K{-in,q) 



and 



B = - 



De^ [ d'^qduodVL 



2tt 



tanh — - coth — - 
{2TT)-i 2T 2T 



{aq - iLU.) 



jK{-in,q) 



(55) 



(56) 



(57) 



(58) 



IV. TRANSPORT NEAR THE QUANTUM 
PHASE TRANSITION 



In the previous section we have derived the expressions 
for the different fluctuation corrections to the conductiv- 
ity that are valid in the vicinity of the entire bound- 
ary (Eq. (|20p ) between the superconducting and nor- 
mal phase in the a-T phase diagram. As the pair- 
breaking strength is increased, superconductivity breaks 
down even at T = once a becomes equal to aco = 
ac(T = 0) = 7rTco/27. In this section we want to map 
out the superconducting fluctuations regimes and obtain 
the finite temperature crossovers that would appear in 
the region 



(63) 



near the pair-breaking quantum phase transition. Note 
however that the validity of our zero temperature results 
will not be restricted to the immediate neighborhood of 



A. Dominant fluctuation corrections 

Since in this section we are interested in the region 
delineated by Eq. (j63p . we will use the asymptotic form 
of the digamma function 



Here 



UJ^ — LU + fi/2 



(59) 



and A and B are introduced for future convenience. The 
first term in each of Sa^/^^ and ScF^t^^ comes from the 
evaluation of the first two density-of-states diagrams, 
while the second term comes from the density-of-states 
diagram that contains an extra Cooperon. 

The "Maki-Thompson" correction can be expressed as 



5u 



MT 



5a 



MT 
sh 



(60) 



where 



D^e^ f d'^qduodVL 
1 



tanh 



LU iK{^i^l,q) 
2T sinh'Jl 



while 



' cth 



'A + iB 



(61) 



(62) 



MT 



with A and B defined above. The contribution 5a 
well as the first term in 5a^lf[ is obtained by evaluating 
the first of the two Maki-Thompson diagrams. The eval- 
uation of the second diagram with an extra Cooperon 
yields the second term in 5a^lj^ which is of the same or- 
der at low temperatures, and of lower order at higher 
temperatures. 



^(2 + 



1 aq — iVli/2 



2tiT 



In 



Uq — iVLi/2 
2^ 



(64) 



in the forthcoming treatment. Let us begin by analyzing 
the Aslamazov-Larkin correction. Using Eq. inside 
the expression for 7 as defined by Eq. (|47p . we obtain a 
reduced form 



7(-if7i,q) 



aq - iVli/2 ac{T) 



(65) 



where the second limit is justified a posteriori by the 
behavior of the integral below. We have 



where 



D^e^ f d'^q dfli 



TTTdaiJ (27r)''sinh'|^ 



(66) 



lmK{'-ini,q) 



rJi/(2a,(T)) 



(67) 

is obtained using Eq. ([32l) . Indeed the integral is domi- 
nated by small values of q and Oi and 5affl^ is critical in 
the vicinity of ac{T). That is also the reason why it was 
legitimate to use the reduced form for K. 

At this point the reader should observe that to ana- 
lytically analyze 5afj^ any further, we need to order the 



10 



energy scales, T and a — ac{T). The frequency integral 
can then be carried out to get the limiting form 



for T <C a — Oic{T) and 



d(27r) 



(68) 
3 (69) 



for T 3> a — ac{T). Observe that the integrand in Eq. 

has exactly the same form as the integrand in Eq. 
(|52p corresponding to the Aslamazov-Larkin correction 
in the temperature vicinity above the classical transition 
at TcQ in the absence of any pair-breaking perturbation. 
Whereas (5t(0,T) and T were the parameters in that in- 
tegrand, here they are 5a{oi,T) and ac(T) instead. 

To be able to access the crossover between the two 
limiting forms, Sa^^^f^ and Jcr^^^, we recast Eq. (|66|) in 
terms of dimensionless variables to obtain 



2D^ 

TTdT^ 



2T 



(70) 



where 

F{v) = 



d'^q dni 



q'Ql 



{2TrY sinh2rJi 



{ij + q^f + n^^ 



(71) 



is a scaling function of the dimensionless parameter 
^ a - ac{T) 



(72) 



In the r] 1 limit, Sa^^^ reduces to Eq. while in 

the ?7 ^ 1 limit it reduces to Eq. The numerical 

evaluation of F(ri) gives the behavior Safj^ over the entire 
range of rj. 

Now let us proceed to analyze the Sa^^f^ part of the 
Aslamazov-Larkin correction. After using the asymptotic 
form of the digamma function given by Eq. ([64|l to get 
a simplified expression, we express <5(T^^ as a sum of two 
terms: 



cth 



d'^qdn^ qlK'\^ini,q) 



coth — ^- 
2T 



sign — 



. ^1 
sign — 



(73) 



At any finite temperature the first (difference) term can 
be shown to be sub-dominant as compared to 5(jf^ while 
for T = it vanishes identically. Thus we need evaluate 
only the second term which is independent of tempera- 
ture. To this end, we carry out a bunch of manipulations 
and do the frequency integration by parts, to obtain a 
neat expression which we denote by 



Trd{d + 2)J {2TTy 



d'^q {Dq^fK^{Q,q) 



(74) 



As for the Maki-Thompson and the density-of-states 
corrections, we find that Sa^f^^ + 5a^f^ is sub-dominant 
as compared to the Aslamazov-Larkin correction Saf^^ . 
On the other hand, the same way as we did in Eq. 
([73)) . we can express ^cr^^'^ and S<7^^f[ as a sum of 
a temperature-dependent term that vanishes at T = 
and a second temperature-independent term. The 
temperature-dependent term is again sub-dominant as 
compared to Sa^^^ and after some manipulation, we find 
the temperature-independent part of A and B (see Eqs. 
([57)1 and ^) to be given by 



2De^ f d'^q Dq^K{0,q) 



and 



TT J d{2iTY 
ADe^ f d'^q {Dq^fK{0,q) 



^° iir J d{d + 2){2Tr)d 
We can re-express Bq as 

i>0 — ^0 



Bo 



2De^ 



d'^q {Dq^YK'\Q,q) 



ind{d + 2)j {2-kY al 



such that we have 



<><^0,cth = ' 



(75) 

(76) 

(77) 
(78) 

(79) 

(80) 
(81) 



to yield 



(5(7, 



AL 



DoS 



Q,cth^^<^0,cth' 

Aq - 3Bn 
2 

^2De^ 



^X^^^T 
^(>'^0,cth 

SBq 
2 

d'^q 



(82) 



{Dq^fK{Q,q) 



TT J d{d + 2){2TTY al 

where Aq, Bq and Bq are all positive quantities defined 
above. The expression for K{0, q) is given by Eq. (pT]) 
which is appropriate for T <C aco- 



B. Fluctuation regimes in the a — T plane 

In the previous subsection we have identified and an- 
alyzed the dominant fiuctuation corrections to the nor- 
mal state conductivity in the vicinity of a pair-breaking 
quantum phase transition out of a superconducting state. 
Based on this analysis, the fiuctuation conductivity is 



Sa{a,T) = SffQ^cthia) + 5a^^ 



(83) 
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FIG. 3: Phase diagram showing the vicinity of the supercon- 
ducting quantum critical point realized via a pair-breaking 
quantum phase transition as illustrated in Fig. [1] The bound- 
ary between the classical and the intermediate regime is given 
by T = a — ac(T). The analytical estimate for the boundary 
between quantum and intermediate regimes is given by Eqs. 
(|94|) . H104p . and (|113|) for the case of a nanowire or doubly- 
connected cylinder, a thin film and a bulk superconductor, 
respectively. For clarity, only the case of d = 1 is shown. 
The quantum regime extends to higher temperatures as the 
effective dimensionality of the system is increased. 



where (Jcro,cth(«) is given by Eq. ([82]) and 8u^^{a,T) is 
given by Eq. At the absolute zero of the tempera- 

ture, the correction to the conductivity is given by 



(5cr(a,0) = (5cro,cth(a) 



(84) 



and it continues to be the dominant correction up to 
a temperature scale T'o(a) at which Saf/I^ {a,T) becomes 
comparable. This regime with a negative fluctuation cor- 
rection which turns out to be almost non-critical, is the 
"quantum regime" with thermal fluctuations playing no 
role whatsoever. In this regime, the Maki-Thompson cor- 
rection turns out to be negative and is half the magnitude 
of the positive Aslamazov-Larkin correction; the negative 
density-of-states correction too is of the same order, and 
the total is negative. The presence of the superconduct- 
ing fluctuations thus impedes the flow of current contrary 
to the naive expectation. It is important to note that the 
quantum regime is outside the scope of any theoretical 
approach that does not take into account the electron 
degrees of freedom in addition to the superconducting 
fluctuations. 

On the other hand, the regime defined by T > a — 
ac{T) is dominated by the Aslamazov-Larkin correction 
Safi^ in its hmiting form <5cr>^;j given by Eq. ([55]) : this 
is the "classical regime" . The "intermediate regime" in 
between the classical and quantum regimes is dominated 
by <5cr^^^, which is the limiting form of the Aslamazov- 
Larkin correction appropriate for T <C a — ac(T), and is 
given by Eq. (j68p . The fluctuation conductivity in both 



the classical and the intermediate regime enhances the 
normal state conductivity due to the additional channel 
of transport via the fluctuating Cooper pairs and is more 
critical than that in the quantum regime. 

The presence of these three regimes means that the 
conductivity behavior depends on how the quantum 
phase transition or the low temperature transition is ap- 
proached. If the quantum critical point is approached 
by coming down in temperature at fixed value of a = 
ttco, then the measurement trajectory sweeps exclusively 
across the classical regime and diverges as the tempera- 
ture tends to zero. On the other hand if aco is approached 
by tuning the pair-breaking parameter at a fixed value of 
T = then the path lies entirely in the quantum regime 
showing the characteristic increase in resistance. On the 
contrary, going away from aco results in a decrease in the 
resistance as a is increased; if the pair-breaking strength 
were to be tuned by a parallel magnetic field, a negative 
magnetoresistance is observed. The intermediate regime 
is swept only while the measurement trajectory crosses 
over from the quantum to the classical regime or vice- 
versa. A trajectory in which the temperature is varied 
at a fixed value of a > aco or a is tuned at a non-zero 
temperature, results in a non-trivial conductivity with a 
non-monotonic behavior. 



1. Nanowire or doubly -connected cylinder (d — 1) 

Consider a pair-breaking quantum phase transition in 
a superconducting nanowire driven by tuning a magnetic 
field, the concentration of magnetic impurities, or yet an- 
other pair-breaking perturbation. As long as the diame- 
ter of the wire is smaller than the superconducting coher- 
ence length, the system is effectively one-dimensional as 
far as superconducting fluctuations are concerned. One 
could also consider a doubly-connected cylinder instead 
of a wire, as discussed in Sec. II. 

For d — 1 the fluctuation conductivity in the classical 
regime (T > a—ac{T)), given by Eq. (|69p . can be written 
in terms of a dimensionless integral to obtain 



De^T 



2-Ka, 



3/2 



dx- 



T 



4\/2 (c 



,(r))3/2 



(85) 



The temperature dependence is more revealing if one uses 
oiciT) ~ ttco — 7r^r^/(6aco) in the vicinity of (see Eq. 
([22]) ). In particular, if the quantum critical point at acO 
is approached by coming down in temperature, then 



6^>,sh{aca,T) = — — 



(86) 



shows a quantum critical divergence with the power T ^ 
when T ^ 0. 



It is instructive to compare the result in Eq. ((85|) in 
the classical regime with the fluctuation conductivity 



<5al^(0,r) = 



cO 



32\/2 [T - T^of'^ 



(87) 



in the temperature vicinity above the classical transition 
at Tco in the absence of any pair-breaking perturbation. 
Derivation of this result requires evaluating Eq. ([52]) 
which contains the same dimensionless integral as in the 
classical regime near the quantum phase transition. 
By using Eq. ([55]) we have 



{a,T) 



5/2 



3a 



dx- 



T2 



12V2 (a-ac(r))5/2 



(88) 



in the intermediate regime between the classical and 
quantum regimes. The fluctuation conductivity contains 
an additional power of T/ {a — ac(T)) as compared to the 
classical regime. 

Expressed in terms on the scaling function that we 
introduced in Eq. ([70|) we have 



■kVt 



with 



ZLl T7-5/2 
96 ' 



?7 < 1 
?7 > 1 



(89) 



(90) 



where 77 is defined in Eq. ([7^ . The full scaling func- 
tion given by Eq. ([7T|) and the negative temperature 
independent term 5cro,ct/i(a), as discussed below, can be 
evaluated numerically. In this way, we are able to get the 
behavior of the fluctuation conductivity as a function of 
temperature and pair-breaking parameter, in the entire 
neighborhood of the superconducting quantum critical 
point. The exact boundary of the quantum regime in 
the a-T phase diagram can then be identified by tracing 
the curve on which the fluctuation conductivity becomes 
zero while changing sign from positive to negative. 

The temperature independent fluctuation correction 
which dominates the quantum regime is given by eval- 
uating Eq. ([82l) . to obtain 



where 



"^^/a^J-oo ha{x)^lnha{x) 



haix) = l + Sa,{a,0)+x^/2 



(92) 



The correction has no critical divergence at a = aco (i.e. 
Sa{a, 0) = 0) and even the first derivative is non-critical. 
The second derivative diverges as (a — aco)"^^^ thereby 
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FIG. 4: Temperature dependence of the fluctuation correction 
to the conductivity (in units of \/De^) in the vicinity of a pair- 
breaking quantum phase transition at a = Oco = O.SSOTco for 
d — 1 (nanowire or doubly-connected cylinder). Plots are 
shown for different values of {a — aco)/Tco- When the super- 
conducting QCP is approached by lowering the temperature 
at a = QcO, the quantum critical divergence of the conductiv- 
ity is given by Eq. (|86p (not included in the figure). The lower 
panel clearly displays the temperatures To{a) (analyticaUy es- 
timated by Eq. (|94p . at which the correction becomes nega- 
tive, there by signalling the entry into the quantum regime at 
a given value of a. 



giving the first non-analytic term in the expansion (ob- 
tained by integrating twice) around aco- We thus obtain 



Sao^cth{a) - Sao.cthictco) 



(93) 



with fli = 0.386 and 61 = — 4\/2/(37r). By using the 
first term in the expansion, one can analytically estimate 
the boundary between the intermediate and the quantum 
regimes to be 



a - 



7/4 



(94) 



TcQ \ acQ 

In Figs. [3] and [5] we display the plots for the fluctua- 
tion conductivity when the vicinity of the pair-breaking 
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FIG. 5: Fluctuation correction to the conductivity (in units 
of ^TOe^) in the vicinity of a pair-breaking quantum phase 
transition at Qco = O.SSQTco for d = \ (nanowire or doubly- 
connected cylinder). Plots show the behavior of 5(7 as one 
sweeps (a — oica) jTca at a given value of T/Tco-The lower panel 
shows the negative correction at T = and a clear upturn 
of the conductivity which is characteristic of the quantum 
regime. This would correspond to a decrease in resistance (a 
negative magnetoresistance, if ol is tuned by a magnetic field) 
as a is increased at T = and a clearly visible non-monotonic 
behavior at low temperatures. 



2. Thin film (d ^ 2) 

Consider a pair-breaking quantum phase transition in 
a superconducting thin film whose thickness is smaller 
than the coherence length which makes the system effec- 
tively two-dimensional as far as superconducting fluctu- 
ations are concerned. The transition can be driven e.g. 
by tuning a pair-breaking perturbation such as a parallel 
magnetic field or the concentration of magnetic impuri- 
ties as discussed in Sec. II. 

The remarkable fact about the fluctuation conductiv- 
ity in two dimensions is that it is an universal quantity, 
independent of the properties of the material under con- 
sideration. The fluctuation conductivity in the classical 
regime (T > a — adT)), given by Eq. ([55]) . can be eval- 
uated for d = 2 to obtain 



dx- 



e2 T 



At: a — adT) 



(95) 



with a^T) ^ Q!co — TT^T^/ (6aco) tti6 vicinity of aco (see 
Eq. (1^ ). If the quantum critical point is approached 
by lowering the temperature at a fixed a = aco, then we 
have 



3e? ckco 



(96) 



which shows a quantum critical divergence T ^ when 

The fluctuation conductivity Eq. (|52p in the tempera- 
ture vicinity above the classical transition at T^q in the 
absence of any pair-breaking perturbation, contains the 
same form of the integrand as in the classical regime near 
the quantum phase transition therefore giving a result 



<5'^f/r(0,T) 



16 T - T, 



cO 



(97) 



quantum phase transition is explored either by sweep- 
ing the temperature or the pair-breaking parameter (see 
the figure captions for details). As discussed above, a 
quantum critical divergence is expected as temperature 
is lowered by sitting at a = Uco and is not shown in 
the figure. If a is increased starting from aco at a fixed 
value of r = 0, the conductivity increases monotonically, 
correspondingly giving a decrease in resistance i.e. a neg- 
ative magnetoresistance if a is tuned by a magnetic field. 
The conductivity shows a non-monotonic behavior as a is 
varied at a fixed T ^ or T is varied at a fixed a > ctco, 
the behavior being more distinct for smaller values of T 
and a — Ofco, respectively. 



analogous to Eq. (1951) ^^'^ consistent with the classic 
literature^. 

By using Eq. (jHS]) the fluctuation conductivity in the 
intermediate regime is given by 

712 2 poo 3 

(98) 



18 ia~adT)y 



It contains an additional power of T/ (a — ac{T)) as com- 
pared to the classical regime, just as we found for d = 1. 
Expressed in terms on the scaling function that we intro- 
duced in Eq. ((70)) we have 
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FIG. 6: Temperature dependence of the fluctuation correction 
to the conductivity (in units of e^) in the vicinity of a pair- 
breaking quantum phase transition at a = Uco = O.SSQTco, for 
the case of a thin fihn {d — 2). Plots are shown for different 
values of (a — aco)/Tco. When the superconducting QCP 
is approached by lowering the temperature at a = aco, the 
quantum critical divergence of the conductivity is given by 
Eq. (|96p (not included in the figure). Note that we have 
plotted the difference 5a{a,T) — 5a{aco,0) which is always 
positive given that Sa{a,0) is most negative at a = Qco (see 
Fig. El). 
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where h = 1.12292. We can henceforth thus consider the 
difference 5aQ,cth{oi) — 5aQ^cth{oico) that has a convergent 
integral to be evaluated numerically. 

Here again as in the d = 1 case, 5<jQ,cth{oi) as well as 
the first derivative has no critical divergence at a = acO- 
From the divergence in the second derivative, we can get 
the first non- analytic term in the expansion around ctco- 
We thus have 

5<jQ^cth{a) - 5<TQ,cth[acQ) (103) 
= [02 (5a(a,0) +b25a{a,Qf'\n5a{a,{)) + ...] 

with a2 0.070 and 62 = l/(27r2). By using the first 
term in the expansion, one can analytically estimate the 
boundary between the intermediate and the quantum 
regimes to be 
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In Figs. [n]and[7|we display the plots for the fluctuation 
conductivity when the vicinity of the pair-breaking quan- 
tum phase transition is explored either by sweeping the 
temperature or the pair-breaking parameter (see figure 
captions for more details). The behavior of the fiuctu- 
ation conductivity is qualitatively similar to the case of 
d — 1. However the critical divergences are weaker and 
the ratio of the fluctuation correction to the normal state 
conductivity is expected to be lower as well. 
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and the full form of F{rj) can be evaluated by doing the 
integrals in Eq. ([7T|) numerically. 

The temperature independent quantum correction 
which dictates the behavior in the quantum regime can 
be obtained by using Eq. ([5^ to have 
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where ha{x) is defined by Eq. ([92]) . For large x the inte- 
grand goes as l/{x In x) and the integral has a very weak 
ultraviolet divergence which can be isolated by evaluat- 
ing the integral analytically by parts (A is the cut-off). 



3. Bulk (d ^ 3) 

Although the pair-breaking quantum phase transition 
out of a superconducting state in low dimensional dis- 
ordered systems attracts more attention, we will now 
consider such a transition in a three-dimensional bulk 
superconductor mainly for the sake of comparison and 
completeness. One could imagine superconductivity at 
low temperatures being destroyed by magnetic impuri- 
ties, for example. The discussion below parallels the pre- 
ceding analysis for d = 1,2. 

The fluctuation conductivity in the classical regime 
(Eq. ([55)) ) is given by 
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FIG. 7; Fluctuation correction to the conductivity (in units of 
e^) in the vicinity of a pair- breaking quantum phase transition 
at Qfco = O.SSgTco for d = 2 (thin films). Plots show the 
behavior of Sa{a, T) — Sa{aco, 0) as one sweeps [a — aco)/Tco 
at a given value of T/Tco.The lower panel shows the correction 
at T = and a clear upturn of the conductivity which is 
characteristic of the quantum regime. This would correspond 
to a decrease in resistance ( negative magnetoresistance) as 
a is increased at T = and a clearly visible non-monotonic 
behavior at low temperatures. 



and at a = arvi reduces to 



(106) 



which is independent of temperature. This is in stark 
contrast to a quantum critical divergence found in the 
case of one and two dimensions on approaching the quan- 
tum phase transition by coming down in temperature. It 
illustrates the fact that fluctuations are stronger and con- 
sequently play a more crucial role in reduced dimensions. 

Same as for a nanowire and a thin film, note that the 
fluctuation conductivity Eq. ([5^ near the classical tran- 
sition at Tco in the absence of any pair-breaking pertur- 
bation is analogous in form to that in the classical regime 
and is given by 
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Although it is critical, the divergence is weaker by one 



power of (T — rco)~^/^as compared to that in two dimen- 
sions, which in turn is weaker than that in one dimension 
by the same power. 

By using Eq. the fluctuation conductivity in the 
intermediate regime is given by 
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It contains an additional power of T / [a — oiciT)) as com- 
pared to the classical regime, just as we found for o? = 1,2. 
Alternatively expressing 5afy^{a, T) in terms on the scal- 
ing function that we introduced in Eq. (|70p we have 



with 



64'/ 

^ „-3/2 
192 ' 



77 < 1 

?7 > 1 



(109) 



(110) 



and the full form of ^(77) can be evaluated by doing the 
integrals in Eq. ()71|1 numerically. 

The temperature independent quantum correction 
which dictates the behavior in the quantum regime can 
be obtained by using Eq. ((82)) to have 
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where ha{x) is defined by Eq. (I92p . The integral 
can again be regulated by considering the difference 
ScQ^cthict) — Sao,cthicxco) as we did in the case of c? = 2, 
and can then be subjected to numerical evaluation. 

By following the same procedure and in one and two 
dimensions we get the expansion 
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+ 63(5a(a, 0)2 + C3(5„(a, 0)5/2 + ...] 

with as = 0.023 and 63 = -0.061 and C3 = 4V2/(157r2). 
Note that the critical divergence shows up for the first 
time in the third derivative to give the first non-analytic 
term in the expansion. By using the first term in the ex- 
pansion, one can make a rough estimate of the boundary 
between the intermediate and the quantum regimes to be 
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The quantum regime is thus expected to extend up to 
higher temperatures as compared to the case oi d = 
1,2. The fluctuation conductivity however has a weaker 
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FIG. 8: Temperature dependence of the fluctuation correction 
to the conductivity (in units of /y/D) in the vicinity of a 
pair-breaking quantum phase transition at OcO = 0.889Tco, for 
the case of a bulk system (d = 3) . Plots are shown for different 
values of (a — aco)/7cO. When the superconducting QCP is 
approached by lowering the temperature at a = aco, there is 
no quantum critical divergence of the conductivity in contrast 
to what is found for d = 1,2 and in fact the conductivity is 
temperature independent (see Eq. pU6p . Note that we have 
plotted the difference 5a{a,T) — 5a{aco,0) which is always 
positive given that 5a{a, 0) is most negative at a = a^o (see 
Fig. El). 



critical divergence as compared to lower dimensions as 
was found above also for the classical and intermediate 
regimes. 

In Figs. [5]and[n]we display the plots for the fluctuation 
conductivity when the vicinity of the pair-breaking quan- 
tum phase transition is explored either by sweeping the 
temperature or the pair- breaking parameter. The quali- 
tative behavior is similar to that discussed for the case of 
one and two dimensions, but the critical divergences are 
clearly much weaker. The biggest difference in the behav- 
ior can be seen if the quantum phase transition at aco is 
approached by coming down in temperature: as shown 
above, there is absolutely no quantum critical divergence 
and the correction is independent of temperature. 



RELATED WORK, EXPERIMENTS AND 
CONCLUSION 



In our work we have studied the superconducting fluc- 
tuation corrections to the normal state conductivity in 
the entire a-T plane, where a parametrizes the strength 
of a pair-breaking perturbation, caused by the presence of 
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FIG. 9: Fluctuation correction to the conductivity (in units 
of I \fiy) in the vicinity of a pair-breaking quantum phase 
transition at QcO = O.SSQTco for d = 3 (bulk superconduc- 
tors). Plots show the behavior of (5cr(a,T) — (5a(acO,0) as 
one sweeps (a — aca)lTca at a given value of T/TcO.The lower 
panel shows the correction at T = and a clear upturn of the 
conductivity which is characteristic of the quantum regime. 
This would correspond to a decrease in resistance ( negative 
magnetoresistance) as a is increased at T = and a clearly 
visible non-monotonic behavior at low temperatures. 



magnetic impurities or a magnetic field, for example. We 
have been particularly interested in mapping out the fluc- 
tuation regimes in the neighborhood of the pair-breaking 
quantum phase transition from superconducting to nor- 
mal state. Our entire analysis has been carried out within 
the framework of temperature diagrammatic perturba- 
tion theory suitable for dirty superconductors and disor- 
dered systems. 

There have been two previous works which have used 
an approach different than ours. Ramazashvili and 
Coleman^° have used an effective action for the pair- 
ing field and evaluated the Aslamazov-Larkin correc- 
tion to the conductivity using a renormalization group 
analysis^^. They have focussed on the quantum phase 
transition driven by magnetic impurities in two and three 
dimensions for weak (BCS) as well as strong coupling 
superconductors and derived the behavior of the con- 
ductivity when the quantum critical point is approached 
by lowering the temperature. The temperature depen- 
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dence matches with what we obtain is two dimensions 
(see Eq. ([M)) ) but in three dimensions our answers do 
not match. While they have considered only the "clas- 
sical renormalization region" , Mineev and Sigrist^^ have 
investigated the entire region around the quantum phase 
transition, using an analogous starting point based on 
the time-dependent Ginzburg-Landau equations (justi- 
fying the use based on arguments by Herbut^^). They 
suggest pressure as a tuning parameter, consider weak 
coupling superconductors in one, two and three dimen- 
sions, and obtain the conductivity corrections in what 
they call, "classical" and "quantum" regimes. Their clas- 
sical regime is identical to ours, while the answers they 
find in their quantum regime correspond to our interme- 
diate regime. 

The above-mentioned approaches based on using the 
time-dependent Ginzburg-Landau equations with a lin- 
ear dissipative time derivative or the corresponding effec- 
tive action are quite powerful and analogous approaches 
have proved extremely useful in the study of quantum 
critical phenomenai^i^iiH However, the Maki-Thompson 
and the density-of-states corrections that require the elec- 
trons to be present in the theory, are outside the scope 
of these approaches and a microscopic calculation, as 
we have carried out here, becomes essential to analyze 
the role of these corrections. Even the zero-temperature 
Aslamazov-Larkin correction that we find is missed in 
the effective approaches. On the other hand, our analy- 
sis is able to identify the regimes where the contributions 
involving the interplay between the fiuctuating Cooper 
pairs and the electrons are subdominant, thereby validat- 
ing the use of the afore mentioned effective approaches 
in those regimes. 

We want to now mention some works which also fol- 
low a microscopic approach, but for physical configura- 
tions different than ours. Beloborodov and Efeto v^^i^^ 
have proposed the negative magnetoresistance observed 
in granular metals in a strong magnetic field and low 
temperatures, to be originating from superconducting 
fluctuation corrections to the conductivity. Galitski and 
Larkin22 have considered two-dimensional superconduc- 
tors in the presence of a perpendicular magnetic field 
and again carried out a microscopic analysis of fluctu- 
ation corrections taking into account all the diagrams; 
they too find zero temperature negative magnetoresis- 
tance. As we have discussed in Rcf. \2§\, negative mag- 
netoresistance is found also for a thin film in a parallel in- 
stead of perpendicular magnetic field. We think it is quite 
a remarkable fact that a negative magnetoresistance at 
zero temperature is a common feature of all these the- 
ories: the Aslamazov-Larkin correction which is always 
positive, the density-of states correction which is always 
negative and the Maki-Thompson correction which has 
no prescribed sign, conspire in all three theories to add 
up into a total negative correction. Although there might 
be different reasons for getting a negative correction and 
a corresponding negative magnetoresistance, it raises the 
question whether this a universal feature of at least a cer- 



tain class of disordered systems in the presence of a mag- 
netic field/pair- breaking perturbation, with fundamental 
reasons at its heart. 

In our work we have considered the corrections to the 
conductivity coming into effect due to the proximity of a 
superconducting state in the low temperature phase dia- 
gram of disordered systems. It is important for our anal- 
ysis that the disordered conductors under consideration 
are assumed to be in the metallic (as against insulating) 
conduction domain. The fact that quantum corrections 
coming from weak localization effects^^ and electron- 
electron interactions (the so-called Altshulcr-Arono\J^ 
corrections) significantly modify the conductivity from 
its Drude-like behavior, even in this domain, does not in- 
validate our analysis. However in comparing our theory 
with experiments, these corrections should be simultane- 
ously taken into account. The weak localization of elec- 
tron waves is a result of an enhanced back-scattering orig- 
inating from the interference between forward and back- 
ward electron trajectories tracing the same path during 
the course of multiple scattering events. It gives rise to 
a negative quantum correction to the conductivity, given 
by the expressioi*i^ 

^L^-^y (l^c(o,,) (114) 

As is evident from the presence of a Cooperon propagator 
Eq. (j26p . the interference effects are strongly diminished 
by the presence of time-reversal symmetry breaking per- 
turbations, thereby yielding a decrease in resistance (a 
negative magnetoresistance, yet again). On the other 
hand, the Altshuler-Aronov correction is determined by 
the diagrams that include only the Diffuson propagators 
and are not sensitive to the magnetic field or other time- 
reversal symmetry breaking perturbations. Thus the in- 
clusion of additional quantum corrections will not affect 
the predicted negative sign of the magnetoresistance (or 
the decrease of resistance with increasing pair-breaking 
strength, in the general case) at low temperatures. 

The experimental effort in superconducting thin films 
has so far been motivated to a large extent by inter- 
est in the so-called superconducting-insulator transition 
(SIT)^, driven by tuning either the disorder (achieved 
by varying the film thickness) or a perpendicular mag- 
netic field. Questions raised by recent experiments about 
the mechanism and interpretation of the transition, has 
revived interest also from the theoretical side. Experi- 
ments on thin films, observing a suppression of super- 
conductivity in the presence of a perpendicular magnetic 
field at low temperatures, have conventionally been in- 
terpreted within the framework of the field-induced dual- 
SIT (acronym we will use to refer to a theory based on 
the boson- vortex duality in a "dirty boson" modeP", al- 
though in the literature such a theory is implicitly im- 
plied whenever the acronym SIT is used) mainly based 
on the negative temperature derivative of the resistance 
above the critical field and finite size scaling analysis of 
the data. Gantmakher et al*^ have made the case that a 
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stringent analysis of many of these data sets might point 
towards inadequacies in such an interpretation. For in- 
stance, the phase interpreted as insulating could very well 
be a metal with quantum corrections; and existence of 
scaling in a limited region might not be sufficient. To 
address these points, they have made measurements on 
NdCeCuO films and found that the microscopic theory 
based on quantum corrections including the supercon- 
ducting fluctuation corrections as obtain by Galitski and 
Larkin^'' could qualitatively describe the main features 
of their experiment including the negative magnetoresis- 
tance (see the discussion in preceding paragraphs). A 
possible crossover between the two interpretations is also 
suggested. Subsequently, Baturina et al.^^ have made 
measurements of magnetic-ficld-dependent resistance of 
ultrathin superconducting TiN films with different de- 
grees of disorder and again concluded that the scaling 
analysis previously regarded as the main evidence of field- 
induced dual-SIT can in fact be observed also for transi- 
tion from a superconductor to a normal metal with quan- 
tum corrections. 

The natural question to ask is, what happens when 
superconductivity at very low temperatures is destroyed 
by applying not a perpendicular but instead a parallel 
magnetic field for which the field-tuned dual-SIT sce- 
nario does not apply? Although there is an interesting 
set of experiments studying the first-order spin param- 
agnetic transition (e.g. Ref. [1^), the case of relevance 
to us is that of a second order transition. Based on their 
measurements on InO films with variable oxygen con- 
tent, Gantmakher et al.'^'' concluded that the behavior 
in the parallel and perpendicular field-tuned case is very 
similar. The non-monotonic magnetoresistance they find 
in both cases, is interesting given that we theoretically 
find a similar behavior in the case of a parallel magnetic 
field, as presented in the previous section. Parendo et 
alj^ have carried out an experiment on ultrathin bismuth 
films to study the thickness-tuned SIT but in the presence 
of a parallel magnetic field. In the immediate vicinity of 
the transition, what they find is a negative magnetoresis- 
tance behavior. Based on the analysis of their data, they 
argue that perhaps its origin could be found in the neg- 
ative fluctuation corrections to the conductivity that we 
find in our calculation. It does seem plausible that the 
conductivity behavior near the superconducting transi- 
tion tuned by varying the thickness at a fixed value of 
parallel magnetic field might be closely related to that 
tuned by varying the parallel magnetic field at a fixed 
value of film thickness. However a definitive calculation 
catering to the former case still remains to be done. As is 
evident from the expression Eq. (jlip . the pair- breaking 
parameter depends not only on the magnetic field but 
also on the thickness of the film. If the only effect of 
changing the thickness were to tune the pair-breaking 
strength, then the two cases would in fact be identical. 
The complication arises because tuning the thickness re- 
sults also in tuning the disorder strength. Tuning the 
magnetic field offers a way of isolating the pair-breaking 



effect. More recently, Aubin et al;^ have done measure- 
ments of field-tuned SIT on NbSi thin films and inter- 
preted the data to conclude that in their experiment, the 
case of perpendicular field is different from that of par- 
allel, based on the presence or absence of a kink in the 
temperature profile of the critical field, respectively. 

It is desirable to have a systematic experimental study 
aimed specifically at exploring the physics of a pair- 
breaking quantum phase transition in superconducting 
films. There have been few works using a parallel field, 
but they have had a limited goal focussing on the SIT, as 
mentioned in the preceding paragraph. Amorphous (non- 
granular) superconducting films that are thinner than ^, 
but not too thin; and are weakly disordered, with as low 
a sheet resistance as possible; would be necessary to as- 
sure that the quantum corrections are small enough. To 
begin with, it will be useful to observe the finite tem- 
perature classical transition and verify the predictions of 
the fluctuation conductivity in its vicinity. By slowly in- 
creasing the pair-breaking strength and lowering the tem- 
perature, one could approach the quantum phase tran- 
sition. Having identified the right films, measurements 
of the temperature and pair-breaking parameter (tuned 
by a parallel field, for example) dependence of the con- 
ductivity, would afford an exciting possibility of discover- 
ing different regimes in the vicinity of the pair-breaking 
quantum phase transition. The increase in normal state 
resistance due to the presence of superconducting ffuctu- 
ations that we find, is in stark contrast to the intuitive 
expectation and is a purely quantum effect. A clear ex- 
perimental signature of such a characteristically quantum 
behavior in the quantum regime, changing over into an 
increase in conductivity in the classical regime, would be 
an important step forward in the study of quantum phase 
transitions and low temperature superconductivity. The 
manifestation to be expected in the plots for conductivity 
as a function of temperature and pair-breaking strength 
can be found in the previous section. The next reveal- 
ing experiment would be to measure the evolution of the 
conductivity behavior with the change of angle made by 
the magnetic field with the film, ranging from parallel to 
perpendicular case. 

Though most experiments on thin films have focussed 
on perpendicular magnetic field and disorder tuned tran- 
sitions, recently Parker et ali^ have performed measure- 
ments on homogeneously disordered ultrathin a-Pb films 
to study the magnetic impurity tuned transition in ad- 
dition. They have compared the conductivity behavior 
near quantum phase transitions tuned by all three mecha- 
nisms. They have concluded that the disorder tuned and 
the magnetic impurity tuned cases show similar behavior 
which seems to be consistent with a fermionic nature of 
the transition to a weakly insulating normal state. Their 
experiment is of relevance to us since they have success- 
fully traced the transition line and shown that it satisfies 
the Abrikosov-Gorkov suppression of Tc given by Eq. ([1]) . 
However the dependence of the conductivity on the mag- 
netic impurity concentration has not been presented and 
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the temperature dependence has been shown only for one 
concentration value that is above the critical concentra- 
tion corresponding to aco- It will be very interesting to 
carry out a detailed comparison of this experiment with 
our theory. 

Systematic experiments exploring the pair-breaking 
quantum phase transition tuned by magnetic impurities 
in three dimensional superconductors would provide a 
good check of our theory, given that in this case we are 
above the upper critical dimension. For example, verify- 
ing the lack of quantum critical divergence as the critical 
impurity concentration is approached by coming down 
in temperature and finding the temperature independent 
behavior instead, would provide a contrast to the critical 
divergence expected in the case of one and two dimen- 
sional superconductors. Finding a weak non-monotonic 
behavior and the presence of a quantum regime would 
be interesting in its own right. In addition, the analy- 
sis of the transition in a relatively simple material would 
provide useful insights in interpreting the behavior near 
the superconducting quantum phase transitions in more 
complex materials. 

As far as superconducting quantum phase transition 
and quantum corrections are concerned, a lot of theo- 
retical as well as experimental work on disordered thin 
films has been carried out over the period of last two 
decades. Relatively less work on similar lines has been 
done on superconducting wires. However the technologi- 
cal advance allowing for the fabrication of superconduct- 
ing nanowires that are uniform and ultrathin (< 10 nm) 
(see e.g. Ref. [11]) has invigorated the field and opened 
up new possibilities. By coating carbon nanotubes or 
DNA molecules suspended over a trench, with a super- 
conducting alloy such as MoGe or Nb, one essentially 
obtains superconducting wires connected on two sides to 
thin film electrodes of the same material. So far most 
experiments done on this set-up have focussed on un- 
raveling the physics of phase-slip fluctuations that result 
in a non-zero resistance below the superconducting tran- 
sition. In one recent experiment, Rogachev, Bollinger 
and Bezryadin^ have looked at the effect of magnetic 
fields. They were able to explain the suppression of the 
transition temperature in terms of the pair-breaking the- 
ory (Eq. (HO])) if the Zeeman pair- breaking effect in the 
presence of spin-orbit coupling was taken into account in 
addition to the orbital effect of the magnetic field applied 
perpendicular to the wire (see Eq. (fT5| V Although the 
fluctuation effects above the transition were not explored 
until now, we believe that the experimental set-up is ap- 
propriately geared to be able to test the predictions of our 
theory. Systematic study of low temperature fluctuations 
both below and above the transition would be essential 
not only to understand the nature of mesoscopic super- 
conductivity but also to be able to successfully control su- 
perconducting electronic circuits involving ultra-narrow 
wires. 

The parallel- magnetic field tuned quantum phase tran- 
sition in doubly-connected superconducting cylinders 



that we discuss in Sec. II -with the pair-breaking pa- 
rameter given by Eq. p^ - has been motivated by its 
experimental realization by Liu et al42,. They clearly 
seem to have observed the enhancement of the conduc- 
tivity resulting from the positive fluctuation contribution 
in the classical regime. But so far there has been no ex- 
perimental evidence of the corrections expected in the 
intermediate and quantum regime. Further experiments 
in this direction would be very interesting. 

In conclusion, we have evaluated the fluctuation cor- 
rections to the electrical conductivity in the vicinity of 
the pair-breaking quantum phase transition using dia- 
grammatic perturbation theory in disordered systems by 
correctly incorporating the quantum fluctuations within 
the formalism. Amongst the three distinct superconduct- 
ing fluctuation regimes that we find, the quantum regime 
is the one is which the contributions to the conductivity 
coming from the interaction between the fluctuating pairs 
and normal electrons are important, while the behavior 
in the classical (or quantum critical) and the interme- 
diate regime is dominated by direct transport via the 
fluctuating pairs. Our theory thus seems to indicate that 
the effective bosonic action formalism should be applica- 
ble outside of the quantum regime. It should be noted 
that even at zero temperature, our theory is expected to 
be valid only outside the quantum Ginzburg region for 
systems below or at the upper critical dimension. One 
important open problem is to understand in what way 
the zero temperature conductivity above the transition 
connects with that below the transition. 

Within our microscopic theory, we have been able to 
extract the finite temperature crossovers to be expected 
near the superconducting quantum critical point. We 
have used the electrical conductivity as a means of prob- 
ing the effects of superconducting fluctuations, however 
a similar analysis could also be carried out for the dia- 
magnetism and other thermodynamic quantities. Exten- 
sion to the case of anisotropic superconductors such as 
high temperature superconductors would also be interest- 
ing. In d-wave superconductors, disorder acts as a pair- 
breaker and the behavior near the corresponding pair- 
breaking quantum phase transition is likely to have some 
resemblance to our findings. 

We believe that pair-breaking quantum phase transi- 
tions form an important class of quantum phase transi- 
tions that not only allow for a systematic, well-controlled 
experimental exploration but also provide the possibil- 
ity of a thorough and comprehensive theoretical analy- 
sis. We hope that our work has amply demonstrated the 
latter and will in turn motivate the former. In the end, 
we expect that the microscopic theory of the supercon- 
ducting quantum critical point and the corresponding ex- 
perimental analysis would serve as a useful prototype for 
understanding quantum phase transitions also in other 
classes of correlated systems. 
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